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Karlqvist Field Approximation 


Head Gap is small compared with gap depth and pole dimensions 
Distribution H,(x,y) is given by circular lines passing through the gap 
comers | 
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Recording Medium: Hard Magnetic Material 


° Hysteresis Loop of the Magnetic Material 
¢ Recording: Field H exceeds Coercivity, after the field is removed, 
- magnetization returs to the remanence level 
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Transitions: Williams-Comstock Model 
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° Magnetic transition is not ideal, it has finite extent( transition parameter (a), 
determined by medium hysteresis squareness (dM/dH), head field gradient 
(dH,,.aq/dx) and medium demagnetization fields (dH,/dx) 


i 
Realistic Transition M(x) 
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Read-Back process: Linear Superposition 


- Head Signal 
Head Gap g (Top view) 


© Head reacts on the changes of medium magnetization 
_¢ Even if the transition is given by ideal step, head integrates the read-back signal (head 
response) | | 
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Read-Back Process 


\ / \ 


Magnetic Medium 


7 
Realistic Transition M(x) 


dH read(x)/dx 


7 an, } 


* Read-back signal is given by the convolution of the medium 
_ magnetization with the derivative of the reading head response (i.e. 
head read-back signal from the ideal step transition) 
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Inter-Symbol Interference (ISI) 


Readback 
signal = 
pl+p2+p3 


read) 9 If (x) = ¥ Area) Sai Mrread ) 5 4. (x) =. p,(x) 
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¢ Read-back signal from adjacent transitions equals to the sum of the 
__ read-back voltages from individual transitions ( linear superposition) 
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Linear ISI: Amplitude Loss 


Dipulse 


«© Dipulse signal: 1 - PW50/T=0.3, 2-PW50/T=1: 3-PWS0/T=2, 4- 
PWS0/T=3; 
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Demagnetization Field 


¢ Static magnetic fields, directed approximately opposite to the direction 
of the magnetization oe 

: Demagnetization fields exist only inside a magnetic material 

¢ They are generated by any discontinuity of magnetization 
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Always opposite to the magnetization 
Infinitely sharp transition: Hd is infinity in the transition center 


Realistic transitions: Hd equals to zero in the transition center 
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Demagnetization fields 


-¢ Far from transition center ( no head imaging) 


M6 
Hj (x) =- 
oo 7X 
with head imaging: 
Ho (x) =-4M5_4 
ff x 


* No head imaging: drops as I/x — 
¢ With head imaging: drops as 1/x? 
¢ Proportional to the medium magnetic moment 
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Hard Transition Shift 


Square wave recorded on DC- 
erased medium 


Easy Transitions are written 
along with medium 
magnetization 

Easy Transitions are not 
distorted 


Hard Transitions are written 


against medium magnetization 


Hard Transitions are delayed 


Hard Transitions may have 
different shape 
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Writing of Hard Transition 


Head 
Movement 


Transition is to be 
written here as 

determined b 
Initial Write 
Bubble (dashed 


Secondary writing 
zone is formed here 


Actual Hard 
Transition is 
written here - at 
the edge of 


smaller write 
bubble 


New write Bubble: --.-"- 
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Writing of Easy Transition 


¢ No demagnetization field, write bubble is not affected 


Write Head 


Gag Head 
Movement 


ST 


Previous Data Magnetization 


eed 


Transition is to be 
written here as 

determined by the 
write bubble 


Actual Easy 
Transition is 
written here - at 


ag the edge of the 
, - write bubble 
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Hard Transition Shift and Magnetic Fields 


Hard Transition : 


Head ee 
cee | Hhead i. written here: H(head)-Hdemag>Hc 


here 
Initial transition 
should be here — 

Hc 
a) 3 


H head - H demag 


H demag 
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Transition Shift 


e Transition is written when H,(x,))=H, 

© Shift: H,(x,+* )=H,+H, 

¢ H,(x,+* )=H,(x))+° dH,/dxl,_ .. = H.+H, 

2 ¢ = -H,/(dH,/dxl_,,) 

¢ Higher density (H,) - larger transition shift 

e Poor writing (e.g. poor gradient) - larger transition shift 

Theoretical equation: - 

, - MO(d + 6/2) 
TQH.r 


M, is remanence medium magnetization, 6 is the medium thickness, d is head-medium 
separation, Hc is the coercive field of the medium, r is the actual size of the head write 
bubble and Q is a numerical factor depending on the head geometry. Function F(d, 6,r) 
depends on head imaging and can be anything from unity to F(d,6,r)=(d+ &/2 ia 


F(d,0,r) 
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Asymmetry Test 


T(p-n) Fw) 


¢ Write Isolated Transitions on DC-erased medium 

¢ Measure distances between peaks positive-negative-positive 
* Asymmetry = Average(T(p-n) - T(n-p)) 
¢ T(p-n)=T+A and T(n-p)=T-A. Therefore: 

* Asymmetry = 2A 
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Overwrite Test 


Old frequency f1 


meas i\ 


New Data 
frequency f2 


Lar eaenenenenee Kiwmmenenme ns 


e A pattern with low frequency f, is written and its average amplitude A1 is 
measured 


_e A pattern with higher frequency f, is written on the same track above the old 
pattern | 
e A residual signal at frequency f, is measured using a bandpass filter or a 
spectrum analyzer. This pattern has amplitude A2 
e The overwrite ratio is calculated as 20log(A2/A1). This value reflects the 
“ability” of a new data to “suppress” old data which is previously written on the 
magnetic medium. Typical overwrite values should be kept below -30 dB. 
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Origin of Overwrite: Hard Transition Shift 


i 
i 


Old 
data 


Head 
movement 


| 
| 
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¢ New pattern with 2f frequency is written on top of If pattern 


e Sequence of 2 easy and 2 hard transitions, - old frequency is 
transformed into the sequence of hard/easy pulses 
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Overwrite ratio =2: Theoretical analysis 


: . ¢ For OW ratio 2 the derivation is 
Assume that the hard transition has shape A(t) and the easy transition has shape e(t). 
Transitions of the old pattern are written with period 47 (the distance between positive an alytic ally stral ghtforward 


and negative transitions is 2T) and transitions of the overwrite pattern are written with 


period 27. All transitions in the old pattern are easy and, therefore, the period of the old ; é ‘ 
signal is given by. | ¢ Overwrite is determined by the 
fold) = et) e(t- 27) 6.5) ratio of hard transition shift A to 
The corresponding period of the overwrite pattern is given by: | the period T of the overwritin 2g 
f (ow) = h(t - A) h(t - T - A) + e(t - 2T) - e(t - 37) (5.6) pattern 


where A is the corresponding hard transition shift. When Fourier transform is calculated, 
each timing shift results in the additional phase of the signal, therefore:. 


F( f(old)} = E(w) - E(w)e"®™" 

F( f(ow)) = H(a)e ~ Haye" + Baye - E@yet | OW 20 | : ei A 
Here w is the frequency which for the old information pattern equals: w;= 27/4T, E( @) is 4 T 
the spectrum of the easy transifion, H(@) is the spectrum of the hard transition. At 
frequency w, complex exponents are having integer number of pi, therefore they result in 


one or i and we obtain that the amplitudes of the old pattern Al and of the overwrite 
_ pattern A2: 


AH) +E) | ¢ Meaning: timing modulation 
A2 =I(1-—-i)( H(@, )e'*" - E(a, ) | OF) index (A IT) 


If we assume that the shapes of the easy and hard transitions are approximately the same, 
the overwrite value equals: 


OW = 2ot04{ icone] (5.9) OW=20 dB 
| _ : ¢ Hard Transition shift should be 
minimized 


¢ Example: T=10 ns, A =Ins, 


If the value of A is small, we can use approximation cos(x) =1 -x?/2, therefore 


Ora) 
T 


OW = 20108 
4 


(5.10) 
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Frequency dependence of Overwrite:Overwrite different LF 
pattern with the same HF pattern. Experiment demonstrates more 
than 10 dB difference between overwrite frequency ratio 10 and 2 
10 20 30 40 50 60 70 80 90 100kfci 


HF pattern density 
(kfc1): 
1-50 
2-66 
3 -83 
4-100 
5-133 
6-150 
7-180 
This is after edge 
trimming, Le. 
effects of track 


width are 
eliminated 


br 
= 


Moat sw oO 
: = 


th 
wo 
=a 


Z om Ue oe ae ee EX AN. ¥ hy v 

OA A oe Of 
t, ‘a Sek ES fa. OB “4 

t_' i} maVilyaw, A ee 
LP ead AN TArAdA 
— AQ ay ohare a ay Att 
A 

| DOK AVANT == a 


_ DENSITY OF OVERWRITTEN PATTERN KFCI 
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Overwrite for different frequency ratios. Theory 


Period of old pattern 


a+d -a+d atd -a+d a+td -a-d a-d -a-d a-d -a-d a+d -at+d a+d 


Easy =a-d Hard =atd 
¢ No simple rule: Sequence of Hard/Easy transitions (e.g. HEHEHHEHEH ...) 


¢ Use ‘average’ (a) and ‘difference’ (d) waveforms of hard and easy 
transitions. Analyze overwrite pattern 
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Overwrite for different frequency ratios. Theory 


Let us use the following trick: imagine that average pulse of voltage a(t) is the average of 
~ easy and hard pulses, including hard transition shift and possible shape difference, and 
_ d(t) is their difference: 


a(t) — A(t) + et) 
a 
d(t)= nN (5.11) 


Using (5.11), h(t)=a(t)+d(t) and e(t)=a(t)-d(t). If we now write the sequence of easy and 
hard pulses in Figure 5.10 using sum and difference of a(t) and d(t) and alternate 
polarities we will see a wonderful picture. Average pulses, as it is expected, behave like 
undistorted sequence - they alternate +a -a +a -a etc. Difference pulses d(t) are exactly 
periodic with the period of the old pattern, indeed from Figure 5.10 we get: 


+d+d+d+d+d-d-d-d-d-d ... etc. 


This is a general rule which holds for any ratio of overwrite frequencies and allows us to 
analyze the overwrite signal. 
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Overwrite for different frequency ratios. Theory 


, the sequence of the average pulses is not modulated by the old pattern and will have 
zero contribution at the frequency of the old pattern. Let the ratio of high frequency to 
low frequency signal be N, i.e. f(old) = 27 /(2NT)=2/ NT. Assume that transitions of 


the low frequency pattern are easy and, therefore, one period of the old signal is given by 


the function: 


f(old) = e(t)- e(t- NT) (5.12) 


As a general rule, the corresponding period of the overwrite pattern is modulated by the 
period of the old pattern and its erroneous part is given by: 


2N-1 


fiow)=¥, d(t—mT)- ¥. d(t—mT) (5.13) 


Only this signal is of interest because alternating average signal will not create harmonic 
at overwrite frequency. When the Fourier transform is calculated at overwrite frequency 


of 2/NT , we get: 
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HF gion | _ ole 2 
F{ f(old)}= co = Ee = ec NT 
- (5.14) 
F{ f(ow)} _ eo ; Sen Par. i PoP 
TNT] I ~ sin(z/2N) 
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Overwrite for different frequency ratios. Theory 


Tf we assume that the shapes of the easy and hard transitions are approximately the same, 
the difference between easy and hard transition is only in the timing shift A and 


| x 
es 
; De) = ose ytexp-i a) - 1 — 1—cos( os A) (5.15) 


the overwrite value equals: 


| 1% 
| 1- cos( A ) 
OW = 20log a +t — * (5.16) 
| sin ( ) | 
2N 


According to (5.16) the value of overwrite should be almost independent on overwrite 
ratio. Indeed, when A is small, cos(x) is approximately 1-x’/2 and, since sin(m/2N) is 
_ decreasing with increasing N almost as 1/2N, we get: 


OW = 20108 =) | | (5.17) 


Theoretical prediction: OW value should not depend strongly on 
the frequency ratio! _ 


April 26 1997 | oe _ A.Taratorin 26 
ST Short Course, SCU 


Overwrite : Proximity Effect 


Gap Movement 


Alread Ree 8 


overwrite pattern Se He 
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Transition 1s to be 
 Iwritten here 

Initial Write 
Bubble (dashed) 


Old transition was written here 
and creates Hd2 which 1s 
opposite to Hd1 


slightly modified 
write bubble 
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When a new 
transition of the HF 
pattern is written in 

the proximity of the 

old transition of the 

LF pattern, the 

demagnetization 

field of the old 
transition reduces 
the effective value 
of HTS. This effect 
is Statistical and 
becomes negligible 
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~Overwrite including Proximity effect 


OW = 2010819 (24 (@ ) 


Where A is the value of HTS, T - period of the HF pattern, ‘¥(@) - 
a function, describing attenuation of HTS caused by transitions of old | 
pattern (decreasing the effective value of HTS at high densities of the LF 


pattern) 


| © ‘¥(@) depends on the density of the LF pattern. This 
function should not depend on the density of the 


overwriting (HF) pattern. 
¢ Proximity effect explains only part of the Overwrite degradation 
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Spectrum of the Overwrite Signal 


P(@) P(o) 


3M) 4@) 5Wy 


° Left: without overwrite, right - when old information is overwritten 


¢ Overwrite causes appearance of side bands around even harmonics of 
the signal at frequencies: 
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Overwrite: References 


Hard Transition shift and overwrite are discussed in: 


N.Bertram Theory of Magnetic Recording, Cambridge University Press, 1994. 


G.Lin, Y.Zhao and H.N. Bertram “Overwrite in Thin Film Disk Recording Systems” - JEEE Trans. 
Magnetics, vol. 29, 6 pp. 4215-4223 - original paper describing relation between overwrite and Hard 
transition shift. : 


_ Proximity effect was studied by Y.S.Tang and X. Tsang in “Non-Linear Transition Shifts in Magnetic 
Recording due to Interpattern Proximity Effects” - Journal of Applied Physics, 74,5, pp.3546-3550, 1993 


The dependence of overwrite on frequency is studied by J.Fitzpatrick and X.Che “The dependence of 
overwrite on proximity Effect” - presented at Intermag-96, Seattle, WA, april 1996, JEEE Trans. Magnetics, 
1996 ( Proceedings of Intermag 96 issue) | 


The spectrum of overwrite signal was obtained by Y.S.Tang and C. Tsang “Theoretical Study of the 
~ Overwrite Spectra due to Hard Transitions effects” - IEEE Trans. Magnetics, vol. 25, 1, pp.698-702, 1989. 
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Non-Linear Transition shift (NLTS) 


Already written 
Area 


>> 
a Region which determines Hard 
Transition shift | 


¢ Origin of NLTS is of the same physical nature as the origin of Hard 
transition shift. The only difference is that NLTS is caused by 
preceeding transitions of the currently written pattern, while HTS is 
determined by the old information. 
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NLTS: Demagnetization field of the previous 
transition 


Write Head 


} 
Movement 


a 
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Previous Transition 
which causes NLTS 


Actual Transition 
is written here - at the 
edge of expanded 
write bubble 


with NLTS=A 


Transition is to be 
written here as 
determined by 
Initial Write 
Bubble (dashed) 


¢ Demagnetization field of the previous transition coincides with the head field direction and expands the write 
_ bubble. Transition is written “early”. NLTS is typically stronger than HTS since the size of write bubble is 
larger than the distance between transitions 
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NLTS: Magnetic fields and equations 


Initial transition 


Head moves 
should be here at Xp H head 


here 
Transition with NLTS =A ; 
is written here: , - 
H(head)+Hdemag>He Tw ee hs 
| ; | 
: A a wie, 
' ome Hiemag + Hh xo+A)> He, 
| A *o 7 
Previous | _ H demag 
Transition 7 
and demagnetization | dH, | 
field | Ay 778 


- NLTS is proportional to the magnitude of the demagnetization field 
_ and inversely proportional to the head field gradient at the transition 
- location 
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The transition is written at the location xpt+A, at which holds: 
Ademag + H, (Xo+A)> He (6. 1) 


where Xo is the location of the ideal transition without NLTS, A is the amount of NLTS. 


Assuming that the value of NLTS is small, we can use Taylor series expansion around Xo 
and obtain that the value of Ais proportional to the value of the demagnetization field: 


demag 


A= 


| (6.2) 
7 


(a 


_ The theoretical expression which allows us to estimate the value of NLTS is similar to the 
equation which was given in Chapter 5 to calculate the value of a Hard transition shift: 


A= aol F(d,0,B) (6.3) 
m0H .B 
where M. is remanence medium magnetization, 6is medium thickness, dis head-medium 
separation, Hc the coercive field of the medium, Q is a numerical factor depending on the 
head geometry and B is the distance between transitions. Function F(d, 6B) depends on 
head imaging and can be unity to F(d,6,B) = (d+6/2)* / B’. The dependence of NLTS 
‘on the distance between transitions is usually described by some power of distance: 
A = N(B) = K/ B™, where the power p is experimentally measured. Typical values of W 
are in the range of 1.5 to 2.5. 
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NLTS: Effect on Read-Back signal 


me 30 0 
TIME.ns 
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I - as expected from linear 
superposition 


2- 25% of NLTS 


Amplitude loss - mainly 
due to stronger ISI 


Second transition is 
moved early 


These distortions are 
catastrophic for PRML | 
channels - samples are 
wrong and Maximum 
Likelihood detector makes 
a lot of errors. 


35 


-NLTS: Write Precompensation 


Correct location X, of the next Transition is 
here, but if written, it will be attracted by 
the previous transition and written earlier 


Previous 


Transition 1s 
written here 


Precompensated Transition 
_|is written with delay p. It 
will: be attracted to the 
previous transition so as to 
become actually written at 


correct location Xo 


e¢ Delay transition during write process if another transition was written 
one clock early 


¢ Total shift of the second transition should be eliminated 
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Precompensation parameter 


NLTS, ns 
By + p— N( Bo + p) = Bo 
p= N(Bo+ Pp) 45° 
6 | a 
pear PA p= N(Bo + p) 

4- ha nin ae ee ST . ge : 
; = 
| , we 

0 : : 

18 20 #22 °&« 24 


Transition Separation, ns 


¢ Optimal Precompensation is smaller than the NLTS value 
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Optimal precompensation and NLTS 


OPTIMAL PRECOMPENSATION 


eSvzoamnd 
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-NLTS and Hard Transition Shift: Hard/Easy 


e First transition 
is written late 
due to HTS 
(A2) 

e Second 
transition is 
easy, shifted 

~ early by NLTS 


¢ Resulting shift: 


April 26 1997 


Dibit 


; 
Movement 
amen at 


een aeeeatt —_—i 


2 @ j#@ ee e# s 


Previous Transition 
which causes NLTS 


Transition is to be 
written here as 


write bubble 
with NLTS=A 


Bubble (dashed) 
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NLTS and HTS: Easy/Hard Dibit 
Write Head 


Secondary writing Zone 
is formed here causing 


Hard transition shift A2 


e First transition 
is easy, not 
affected 

¢ Second 

‘transition is 
pulled early by Nee 
NLTS and late _ 
by HTS 
¢ Total shift: 


= @ mm @ eae 8  «@ 


Previous Transition 
which causes NLTS 
Al was Easy 


Actual Transition 
is written here - at the 
edge of modified write 
bubble 


with NLTS=Al-A2 


Transition is to be 
written here as 
determined by 
Initial Write 
Bubble (dashed) 


Anja =A, —Ap 
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NLTS and HTS: Easy/Hard Dibit 


The resulting write bubble is determined by the balance of two demagnetization fields 
and the actual transition will be written at the location, where: 


H(head) +H1(demag)-H2(demag) > Hc (6.6) 


Similar to the analysis presented in section 6.1, if both demagnetization fields are small 
enough, the resulting shift of the transition: 


Al ~ 42, 
A eet lemag lemag | s 
(# | nis 
dx ** 
or, equivalently: 
A rin =A,—A, (6.8) 


where A, is the value of NLTS without the hard transition shift and A, is the value of 
the hard transition shit. Therefore, for an Easy/Hard (E/H) dibit transition, NLTS and the 
hard transition shift have opposite directions and resulting transition shift will be smaller 
than for AC-erased medium. 
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Hard/Easy and Easy/Hard Dibits: Summary 


| Apu 
| ae | ; AgnvzAcA 
F/H dibit | | I+ ai Aran 


Figure 6.10 Interactions of NLTS with Hard Transition Shift 


April 261997 | | A.Taratorin 
| LIST Short Course, SCU 


42 


_ Hard/Easy and Easy/Hard Dibits: Read-back 


NEGATIVE DC POSITIVE DC | 
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Amplitude Loss ( Partial Erasure) 


Transition with 
Partial Erasure 


. Isolated | | ~ Readback after 
pulse | partial Erasure: Isolated 
| | pulse 


| Linear Read-back 
signal is dashed 


¢ Transition percolation causes amplitude 
loss 


co oe = -Read-back signal: Veg (x) = (1— V(x) 
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Spectral Measurement of Partial Erasure 


e Assumption: while PE and NLTS coexist at high density, if a square 
wave pattern is written, all transitions are shifted equally and, 
therefore, NLTS does not interfere with measurements 


- Original transitions 


| 
pT : 
<—> 7 1 1 
T 
Written transitions 
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Spectral measurement of Partial Erasure 


If a periodic pattern consisting of pulses of alternating polarity with distance between 
pulses T is written, the spectrum of this pulse sequence consists of odd harmonics at a 
fundamental frequency of@, = 272 /2T. The shape of this spectrum repeats the shape 
P(q@) of continuous spectrum of an isolated transition p(t). The Fourier transform of this 
readback signal is given by: | 


Ryp(@) = @ P(@) ¥ 5(@ — ka) (6.11) 


k=odd 


Note that the spectrum amplitude is proportional to the pattern frequency @ =2/T. Let 
us write two patterns: low frequency (LF) and high frequency (HF) with different 
distances between transitions. Now it is easy to demonstrate that in general, at any odd 
harmonic k ( periodic pattern does not contain even harmonics) we have: 


| I Myr 
R,-(k@,;) = Rup(@yp), Where O,; = i (6.12) 
Let k=3. The first harmonic of this pattern equals: 
Rup (@)) = O@)P(@) (6.13) 
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Spectral Measurement of Partial Erasure 


Let us now compare what happens if we write a pattern with a lower frequency, for 
example @,/3. In this case we will look at the third harmonic of the pattern which 


equals exactly @,: 
M 
R,-(@)) = a P(Q) (6.14) 


As we see, R,-(@)) = Ryp(@))/3. Since this pattern is low frequency, it will not contain 
NLTS and partial erasure. Therefore, if we measure 3R,,(@,) and compare this value to 
Rij (@)), we can easily estimate amplitude loss: _ 


3R, (Wy (1 — 202) = Rye (@p) (6.15) 


A factor of 20 appears, Since in a high frequency pattern each transition has partial 
erasure from its two neighbors: one to the right and one to the left. Therefore, 


L Ry (@y) 
C= 6.16 
2 6R,-(@) a 
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Partial Erasure measurement: Summary 


e Define LF pattern as 100100100... in NRZI encoding, i.e. each “1” is a transition and 
each “0” is no transition. : 
e Define NRZI HF pattern as 11111.... 
e Perform the overwrite test for the frequency Wo of the HF pattern (Refer to Section 
pyre) Me Moree 
= 1. Write the LF pattern and measure the amplitude V(LF) on the output of the 
band-pass filter at @p (or use a spectrum analyzer) 
2. Write the HF pattern and measure its amplitude V(HF) at @o. 
3. If steps 1,2 are done using a standard overwrite test, the overwrite ratio OW is 
calculated as OW=20log(V(HF)/V(LF)). 


¢ Calculate amplitude loss as 
| 1 V(HF) 
C2 6.17 
2 OV(LF) ( ) 


OF, if OW value in dB 1s obtained, 
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Partial Erasure: Density Dependence 


PARTIAL ERASURE 


35 


30 


RO 
gi 


RO 
Low ] 


Partial Erasure, % 


30 50. 60 #70 #4280 #490 #4100 «+140 120 
KFCI 
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NLTS & Amplitude Loss: References 


¢ A. Armstrong, H.Bertram & J.K. Wolf "Nonlinear Effects in High Density Magnetic Recording" - — 
IEEE Trans. Magnetics, 27, 5, pp. 4366-4376, 1992 - A useful paper describing non-linearities and 
their models. | 


: H.Bertram, A. Armstrong & J.K. Wolf "Theory of Nonlinearities and Pulse 
° Asymmetry in High Density Magnetic Recording" - IEEE Trans. Magnetics, 
© 28, 5, pp. 2701-2706, 1992 - model calculations of NLTS. 


¢ Interactions of NLTS with hard transition shift were first mentioned by C. Tsang, Y.S. Tang " Time 
_ Domain Study of iia Effect Induced menses Shifts" - IEEE Trans. Magn. 27, p. 795-802, 
1991. 


¢ The most comprehensive paper describing precompensation, partial erasure and spectral 
measurements of partial erasure is by X.Che “Nonlinearity Measurements and Write 
Precompensation Studies fora PRML Recording Channel” - IEEE Trans. Magnetics, vol. 31, 6, 
~ pp.3021-3026, 1995. 
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NLTS Measurements 


e¢ Spectral Elimination (5-th Harmonic) : write simple repetitive pattern, 
containing dibits and isolated transitions, measure specific harmonics 
using spectrum analyzer 

¢ Pseudo-Random Sequences (Dipulse Extraction and Time Correlation) - 
write a special pseudorandom pattern, digitize it with high resolution and 
perform mathematical analysis ( Fourier transform or correlation) 
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Harmonic Elimination 


¢ Write a special pattern on the disk and measure the amplitude of a 
particular ( e.g. 5-th ) harmonic of this pattern E(kw,). The pattern is 
constructed in such a way that if the superposition of all pulses is 

linear, the corresponding harmonic is equal to zero. 


¢ Write a different pattern consisting of transitions with a specified 
period on the disk and measure the amplitude of the same harmonics of 
this pattern X(kw»). 


¢ Compute a value of NLTS, usually as 
5= E{kay)/ [key X(key) 
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Principle of Harmonic Elimination 
Periodic pattern y(t), p(t) - isolated pulse 


y(t) = p(t)— p(t-T,) + p(t-T)-..-(-1)" p(t-Ty ) 
- Fourier Transform: 

¥(@) = P(@)|1 =e 1°% +e FOR, -(-1)% eT | 
At k-th harmonic Y(kw,)=0 


 k@) =k2a/ NT 
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Example of Pattern 


N=2m+2n+6 


..m0Q’s... 1...n0’s... 1 1...m0’s...1...n0’s.... 


© 1100..(m 0's)100..(n 0's)1100..(m 0's)100(n 0's) 
¢ A full period of this pattern equals 2m+2n+6. 


i e 
e nprie: M—ph=—Oo QCTIOC IN a 
a a 
7, LZ ~, ~~, , ~ ~~ J SH 7 XS 
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Error Signal from NLTS 


m@ Error signal from dibit with NLTS: 


Pied wl ied 


dt 


Spectrum of Error signal: 


-0.4 


EK ay) = Akon PK) 


Reference signal: 


R(keyy) = P(Kea )(1—exp{—i})= 2P(Kay) 


1 E(k@))/2 1s E(k@p) 
kM) R(kK@_))/2 = =k@ R(kQ,) 
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Pairs of pulses eliminating k-th harmonic 


pr 


— Y¥(@)= P(a)\i+e'” ) 


kp 


-  —120— 
Y(k@))= P(k@))A+e ) 
21 
OQ = -— 
NT 
2m+1 
_ kp = N 
¢  This.pair of pulses will result in zero spectrum at: 2 
|e Wecan combine positive and negative pairs of pulses to create 
harmonic elimination patterns > 
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Simple Pattern for NLTS measurements 


Pattern period N=30T 


° This pattern conisist of: pair of positive pulses (P1 and P3, distance 9T) and pair of negative pulses 
_ (P2 and P4, distance 15T). P1 and P2 form dibit and pulse P2 is distorted by NLTS. Pattern period is 
30T | 


lO ee 
k@) R(K@))/2 kW R(k@) 
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Precompensation Mode 


Pattern period N=30T | PRECOMPENSATION PLOT 


fe 


PRECOMPENSATION, ns 


. Delay pulse P2 and Measure amplitude of 5-th harmonic for different delays. 
¢ Minimum point - optimal precompensation for dibit 
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Bit Period Mode 


Pattern period N=30T Z BIT PERIOD PLOT 
F le 
oT 
PL 7. PS 6Bh- 
2° 
p 4} 
a | 
= 4\.. 
2 _. 
P2 ; 
; ; P4 3 ; 3 3 
S$ $$ >: 3 0 7 7 | 3 
amen 10 12 14 16 18 20 
re | BIT PERIOD, ns 


¢ Delay pulses P2 and P4 by the same amount 

¢ Measure 5-th Harmonic 

¢ Effectively change bit distance between P1 and P2 while keeping elimination property between P2 
and P4 


eet 
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Partial Erasure and 5-th Harmonic measurement 


Let us consider what happens with the read-back signal for a pattern shown in Figure 7.4 


when it contains amplitude loss o due to partial erasure. The dipulse signal distorted by 
_ both NLTS and PE is given by: 


f(t)=(-a@)p(t)-CU-a@)p(t-T +A) (7.15) 


The error signal (the difference between the ideal dipulse and the signal given by (7.15) ) 
generated by this dipulse is: 


d 
e(t) = ap(t)— ap(t-T) +(1- aA (7.16) 
| The spectrum of the error signal is, correspondingly: | 
E(w) = P(w)(a — ae” + iA(—a)@) (7.17) 


At k-th harmonic E(k@,) = P(k@,)(a - ae**"" + iA(- a) kan, ). 
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Partial Erasure and 5-th Harmonic measurement 


Eq. (7.17) is used to calculate plots shown in Figure 7.8. This figure is a dependence of | 
the measured value of NLTS using 5-th harmonic on actual value of NLTS versus 
different amounts of partial erasure (0% to 50% with 5% steps) 


As seen, both NLTS and partial erasure will affect spectral elimination results. For 
example if 20% NLTS and 20% amplitude loss are present in the signal, the spectral 
method will measure approximately 35% of NLTS. On the other hand, if NLTS is absent, 
the measured signal will be completely due to partial erasure. 


_. When dibit separation 1s proportional to the total length of the harmonic elimination 
pattern, k=5, N=30, we can obtain a simplified expression for (7.17), where: 
E(K@,) = P(ka, | — oe 3 + iA(1— a@)ka@, (7.18) 


To find an absolute value we find the squares of real and imaginary parts in (7.18) 
substituting cos(z/3)=1/2 and sin(7/ 3) = ./3 /2. If we discard all cross-terms having 


products of either A’ or a”, assuming that their values are small, (7.18) approximately 


equals: 
7 2 a” 2 
|E(key )i=1PCkeog fo +A + [Aa (7.19) 
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Influence of Partial Erasure on 5-th harmonic 


MEASURED NLTS % 
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Separating NLTS and Partial Erasure 


AA’ +BA+C=0 


2 
z =| Stew) 


P(K@p ) 
A=(1-@)*(k@)” ; 
B =(2(1— a)k@y sin(k@pd) hee Bt vVB° -4AC 
~C=a(1—cos(k@pd))” +a” sin*(k@yd)—Z 2A 


Let d=T, k=5, “0 =2n/30T. The measured value of a=0.2. Amplitude of the error signal 
E(k®)) from a single dibit equals 6mV. Amplitude of the reference pattern consisting of 

& UNS 
isolated pulses with distance between pulses 15T equals 40mV, therefore P(K@,)_90 
mV. We find: . 


Z=0.098, A=0.7,B=0.29,C= -0.058 


From here the estimated value of A=14.8%. At the same time, 
direct measurement will predict almost 29% of NLTS. 
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Influence of Noise on 5-th harmonic 


Spectral NLTS measurement for different SNR 
0.6 


— 605 


. 
. 
Peer eTT ee TTT Tee SETTLE TE Te ee ree ee rrr rr rrr irr rir en de eee mennamennesencdosossgitannces 


SNR=i0 db 


Measured NLTS 


Be a 
‘ _ : i SNR=40 db: : : 


1 
PSST TTT ee Le ee eee eee ee EOL LLL CL 


0 0.0 0.1 015 02 025 0.3 0.35 


0.4 0.45 0.5 
Ratio of NLTS to Bit Period 
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Pseudo Random Sequences 


¢ Write a pseudo-random sequence on disk 


° Sample a full period of read-back signal with high resolution analog- 
to-digital converter (ADC) 


© Process resulting digital signal using either Discrete Fourier Transform 
or Correlation analysis 


¢ Polynomial x’+x3+1: 127 bits sequence y(i+7)=XOR {y(i),y(i+3)} 


10000001 00010011 00010111 01011011 00000110 01101010 
01110011 11011010 00010101 01111101 00101000 11011100 
01111111 00001110 11110010 1100100 
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Shift and Add Property 


dae Product of 2 elements of PRBS is the same sequence with some shift XOR(y(k), y(k+p)) = y(k+M) 
¢ ~~ XOR is equivalent to product of NRZ bits 


m When NLTS is present, it creates the terms proportional to product of 
several bits , e.g. a(k-1)a(k)a(k+1). 
}m > Read-Back voltage is approximately V,(t)+kV, (t+M) 
m M=25.5 ; kis proportional to NLTS 
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Theory of the dipulse extraction method 


Let us write a readback voltage from the head as the following sum: 
V(t) => Salk +1)— aka kD) (7.31) 
; | 


Here h(t+kT) describes the shape of an isolated pulse. Indeed, if pattern is defined in NRZ 
terms, the transition corresponds to a change from +1 to -1 magnetization (negative read- 
back pulse) or to a change from -1 to +1 magnetization (positive read-back pulse). If 
NLTS equal to A is present, readback can be written in the following way: 


V(t)= >) [a(k +1)—a(k)|h(t + kT - lack + 1) —a(k)][a(k) — a(k —1)) (7.32) 


If at least one of the differences a(k+J1)-a(k) or a(k)-a(k-1) is equal to zero, (7.32) 
becomes equal to linear signal (7.31). However, if the current transition is preceded by 
the previous one, the term A is added to the delay of the h(t+kT) and NLTS is present. 
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Theory of the dipulse extraction method 


A similar trick is used to include different non-linearities into Eq.(7.32). For example, 
Hard/Easy transition shift equal to value e€ is described by the term 


—5 + Fla(k +1)—a(k)] For a particular polarity of the transition, the difference of 


a(k+1)-a(k) = 2 and the resulting shift is zero. For a different polarity of transition the 
shift equals to €. | 


The next step in analyzing Eq.(7.32) is to expand the function A(t) into the Taylor series 
- around t+kT . The Taylor Series expansion results in quadratic terms of pseudorandom 
sequence coefficients. After a sequence of mathematical transformations which could be 
found in the original paper by Palmer et.al., it appears that the main contribution into the 
playback signal is introduced by the cross-terms of the pseudorandom sequence. For 
example, the dominant contribution for NLTS is given by the term a(k-1)a(k)a(k+1) . 


For the polynomial x/+x3+1 this term causes the shift of M=-25 bit periods. The 
resulting waveform is i tice written as: 


A A 3T. A A 
MO Vir S) tse vert 7) ap VETS +[M- 310) (7.33) 


As seen Slee Eq.(7. 33), the read-back signal consists of its several delayed and shifted 
versions (echoes). Several echoes are created by NLTS, (main, delayed by 1.5T and 25.5 
bit ease 


April 26 1997 a _ A.Taratorin 


JIST Short Course, SCU 


68 


Dipulse Extraction Method 


¢ Write Pseudorandom Sequence on Disk 


e¢ Digitize and extract read-back signal from exactly one period of this 
sequence | 


¢ Resample Data to obtain N samples per bit period 

¢ Create oversampled NRZ sequence of the same length 

¢ Calculate Fourier Transform of Read-back data 

< _ Divide it by Fourier Transform of the NRZ sequence 

¢ Obtain Dipulse Plot. If NLTS=0, echo at 25.5 bit periods is absent 
= -NLTS/T = 2(amplitude of echo)/(amplitude of main dipulse) 
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Pseudo-random sequence Read-Back 


200 400 600 800 1000 


¢ Exactly one full period of the pseudorandom sequence should be | 


captured — 
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Dipulse Extraction Plot (WITE software) 


Bit Periods 
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Meaning of Dipulse Plot 


Methods for Characterization of Linear Systems 


HEAD SIGNAL 
WRITE DATA Step Response 
| ISOLATED PULSE 
Linear System 
WRITE DATA 
DERIVATIVE OF 
Impulse Response 7 PULSE 
Linear System | 
DIPULSE 


WRITE DATA Dibit Response 


IL 


Linear System 
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Method of Time Correlation 


Correlation Function — 
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Algorithm of Time Correlation measurement 


“a 


¢ Write a pseudorandom sequence 

¢ Digitize at least two full periods of this sequence V(t) and V(t+NT) 
¢ Calculate their cross-correlation function R(u)=R(V(t), V(tt+NT+u)) 
© SNR =R(0)/(1-R(0)) 

° Peak of R(d) at d=25.5T corresponds to NLTS 

* NLTS/T=-2R(d)(1+1/SNR) _ 
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Time Correlation Plot 


10 20 30 40 90 60 
BIT PERIODS 
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Partial Erasure and Dipulse Extraction 


DIPULSE PE=20% NLTS=30% 


24 25 26 
BIT CELLS 


¢ Echo of PE is 1/2T away from NLTS echo (left figure) 
¢ Superposition of echoes is almost linear:1 - PE=20% NLTS= 0%; 2- PE=0%, NLTS=30%,; 3 - 
PE=20%, NLTS=30%; 4 - sum of echoes 1 and 2 (night figure) 


a a a 
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Partial Erasure and Time Correlation 


CORRELATION PLOT 


5 20 
BIT PERIOD 


¢ Partial Erasure Echo is 1/2T away from NLTS echo (left figure) 


¢ Superposition of NLTS and PE echoes is non-linear (right figure): 1- 30% PE and 0% NLTS ; 2- 
0% PE and 30% NLTS ; 3- 30% NLTS and 30% PE 
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~ Mix of PE and NLTS in Time Correlation 
method 


TIME CORRELATION PERFORMANCE 


INLTS=50% 


. 
ry 
’ 
ry 
v7 


MEASURED TOTAL DISTORTION % 


0 5 10 15 200 «26 30 35 40 45 50 
| PARTIAL ERASURE % 
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Influence of PW50 on Pseudorandom methods 


RATIO MEASURED NLTS/REAL NLTS . 


1.5 2 
RATIO PW50/T 


Dipulse extraction is more critical to PW50 than Time Correlation | 
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NLTS Measurements: Summary 


e¢ Harmonic elimination methods are relatively simple to implement 
¢ Harmonic elimination method measures NLTS for a dibit transition 


¢ Harmonic elimination methods may predict optimal precompensation 
parameter and dependence of NLTS on bit period 


¢ Pseudorandom methods are complicated 
¢ Pseudorandom method measure NLTS for random pattern 
¢ AlUNLTS measurement methods are affected by partial erasure 


e If Partial Erasure is measured, it can be separated from 5-th harmonic 
‘Measurement | 
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NLTS Measurements: References 


e Harmonic Elimination methods are described by: 
e YY. Tang, C. Tsang " A Technique for Measuring Non Linear Bit shift" - JEEE 


° Trans. Magn., 27, 6, pp. 5316-5318, 1991. First paper describing method of 

oc harmonic elimination for NLTS measurement. 

e  X.Che, M.J.Peek and J.Fitzpatric "A Generalized Frequency Domain 

o Nonlinearity Measurement Method" - IEEE Trans. Magn. Vol.30,6, P. 4236,1994 
° | - Description of the second method of harmonic elimination 


° X.Che “Nonlinearity Measurements and Write Precompensation Studies for a PRML Recording Channel” - 
_ IEEE Trans. Magnetics, vol.31, 6, pp.3021-3026, 1995 - a description of the harmonic elimination method 
and the influence of partial erasure on harmonic elimination. 


e  Pseudo-random algorithms and dipulse extraction method are discussed in 


- OD. Palmer, P. Ziperovich et al "Identification of Non Linear Write Effects 
Pe 5 Using Pseudorandom Sequences" - JEEE Trans. Magn, 23,5, pp. 2377- 


© 2379, 1987 -Original paper describing Dipulse Extraction method. All mathematics) 1S 
concentrated in the Appendix and is given almost without explanations. Requires significant 
effort to read and follow the authors. 
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NLTS Measurements: References 


P. Newby and R. Wood "The effects of Nonlinear Distortion on Class I'V 
Partial Response’ - IEEE Trans. Magnetics, 22,5, pp.1203-1205, 1986 - First 
paper describing application of pseudo-random sequences to study influence of 
| NLTS on the quality of PRML signal. 
_ X.Che and P. Ziperovich "A Time Correlation Method of Calculating Nonlinearities 
Utilizing Pseudorandom sequences" - JEEE Trans. Magnetics, vol.30,6, p.4239, 1994 


-This paper presents a description of the time domain correlation method, in aaa! 


the dependence of results on PW5O of the pulse is studied. | 

_ G. Mian and T. Howell "Determining a Signal to Noise Ratio for an Arbitrary Data 

Sequence by a Time Domain Analysis" - IEEE Trans. Magn., vol. 29,6, p.3999, 
1993 - Original paper describing the principles of correlation analysis for SNR 

measurement | 

~ G.Mian "An Algorithm for a Real Time Measurement of Nonlinear Transition 

Shift by a Time Domain Correlation Analysis" - IEEE Trans. Magnetics, January 

1995 - Original paper describing the principle of Time domain correlation NLTS 
measurement 


April 261997 a | | A. Taratorin 
a | | | IIST Short Course, SCU 


82 


Interactions between Transitions 


¢ NLTS(2) is large, NLTS(3) is smaller, NLTS(4) is again larger 
¢ NLTS depends on several previous transitions 
e Different NLTS levels for each transition in a series??? 
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Finding NLITS for a Series of Transitions 


N(B) 


AS 


Point of 
convergence Ao 


A2=N(Bo+Al) 
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NLTS settling for a series of transitions 


NLTS Settling 


35 


—  & 


NLTS, ns 


” 4 6 8 10 12 14 
Transition Number 


No 


| k 
A(k) =} (-1)""N[mBy + A(k = m)] 


m=1 


A.Taratorin 
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Dibit NLTS vs NLTS for Random Pattern 


x7+x3+1 polynomial contains total of: 


° 8 dibit transitions | 

a | | 4 tribit transitions ( 3 adjacent transitions) 
° : | 2 series of 4 adjacent transitions 

° | 1 sequence of 5 adjacent transitions 
° - 1 sequence of 7 adjacent transitions. 
° Total we have: 

. 16 transitions with value of (2) 

. 8 transitions with value of (3) 

. _ 4 transitions with value of (4) 

° _ 2 transitions with value of (5) 

° 1 transitions with value of (6) 

* 1 transitions with value of (7) 

° When echo is measured - 


_ 16A(2)+ 8403) + 4A(4) + 2A(5) + A(6) + AC7) 


A 
32 
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5-th harmonic NLTS vs Pseudorandom NLTS 


NLTS PSEUDORANDOM % 
rm) 
a 


~h 
nn 


eeerrere reer rn Serre rere etre etter eet ererereteeeeree sereeteeerls Otte Phere ret erttertrtens fetter rent elon See sent iren enn Sennen eee rene Seen Ceeereree rere 


O° 5 10 15 20 25 30 35 40 45 50 
NLTS FOR DIBIT % 
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Interactions with Hard Transition Shift 


e Easy Transition is affected by previous Transitions 


¢ Hard Transition is affected by previous transitions and by hard 
transition shift | 


¢ NLTS for a series of transitions will never settle 


¢ Transition shift oscillations are larger than the values of hard transition 
shift (T=l6ns, Hard shift= Ins, NLTS=4ns, oscillations are 2.5 ns) 


¢ Oscillations are different depending on whether the first transition in a 
train of transitions is hard or easy 


¢ Serious problem: precompensation is not effective 


Cee ne en ce cs 
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Spectral Measurement of Individual Transition 


Shift 
| \_S\ 
Ve VY 


Only NLTS of the 3-d Transition will be measured 
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Model and Experiment:NLTS versus Transition 
Number 


APPROAIMATION OF NLTS 


15 20 25 30 35 
TRANSITION SEPARATION,NS 


e - Approximation of N(B) dependence used for model calculation 


° O - measured x-model 
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Model and Experiment: NLTS versus Transition 
oe Number 


TRANSITION SHIFT, AC-ERASE 


TRANSITION SHIFT, NS 


A , 6 
TRANSITION NUMBER 
Dependence of NLTS on transition number. 


Medium is AC-erased. o-measured data, 
-X- experimental data 
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Model and Experiment: NLTS versus Transition 
Number 


_. TRANSITION SHIFT, NEGATIVE DC-ERASE 
TRANSITION SHIFT, POSITIVE DC-ERASE | , 


M6 
i 5 ma 
ye w 

z= 
84 S 4 
= a 
2 a 
So r 2 


4 6 8 0 
| — TRANSITION NUMBER 
_ Dependence of NLTS on transition number. 
Medium is positive DC-erased. Dependence of NLTS on transition number. 
O- measured data, x- experimental data Medium 1s negative DC erased. 
| - - O measured data, x- experimental data 


4 6 
TRANSITION NUMBER 
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Model and Experiment: NLTS versus Transition 
Number 


TRANSITION SHIFT, NEGATIVE DC-ERASE 
8 a : 


MODEL, GROUP OF 20 PULSES 


_o 


Bs n6 
~ -~ 
rare rae 
7 am 
she w 
Zz 
5 4 & 4 
A = 
2 D3} - 
Oo »~ = 
F 2 F 2 


exnds 


9) 


: 8 
5 10 
TRANSITION NUMBER | TRANSITION NUMBER 


Dependence of NLTS on transition number. 
Medium is negative DC-erased. 
o- measured data, x- experimental data 
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_ What happens if HTS interacts with NLTS? 


¢ HTS modulates transition shift in a square wave pattern. 
¢ Transition shift oscillates between largerand smaller values. 
¢ Resulting oscillations exceed the value of HTS 


A.Taratorin, J. Fitzpatrick, S. Wang and B.Wilson, Non-Linear Interactions in 
a Series of Transitions, IEEE Trans. Mag, January 1997 


¢ Excessive oscillations of transition shift 

¢ Precompensation is not effective: we don’t know the old information 
pattern 7 
¢ Old story: Overwrite 

_*- Non-linearities limit performance of magnetic recording channels 
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Precompensation of Adjacent Transitions 


Pe = > (-D"™" Nm + py) 


m=1 


¢ Individual Precompensation levels are required for each transition in 
- series 
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Optimal Precompensation vs Transition Number 


OPTIMAL PRECOMPENSATION 


PRECOMPENSATION, ns 


2 3 4 5 6 ¢ 8 
| Transition number 
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Transition Modulation and Head/Media 
Parameters 


NB) 
e¢ N(B) is steep (curve 2): influence of 


previous transitions is smaller, but 
influence of hard transition shift is 
! larger 
, ! 7 ¢ N(B) is flat (curve 1): influence of 
! | previous transitions is larger, 
individual precompensation is 
required but influence of hard 
transition shift is smaller 


eee wlio me me me me ea ee ee oe OE Oe Ne ee ee 
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_ o —— of Overwrite after Edge Trimming 
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7 ‘Overwrite bh isa function of old pattern density 
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Overwrite including Proximity Effect 


A 
OW = 20logig|7-¥ (@) 


Where A is the value of HTS, T - period of the HF pattern, ‘¥(o) - 
a function, describing attenuation of HTS caused by transitions of old 


pattern (decreasing the effective value of HTS at high densities of the LF 
pattern) a | 


e ‘Y¥(@) depends on the density of the LF pattern. This 
function should not depend on the density of the 
overwriting (HF) pattern. | 
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Equivalent HTS - product A¥(a) 


HF pattern density 
(kfc1): 
1-50 
2-66 


3-83 
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6-150 
7-180 


10 20 30 40 50 60 70 80 90 
| DENSITY OF OVERWRITTEN PATTERN KFCI 


Curves 1,2,3,4 (up to 100 kfci) - same trend - good agreement with theory 
| Starting from 100 kfci - amplification of HTS, deviation from the theoretical 
model 
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NLTS and amplification of HTS 


Precompensation fixed at 0.00% of bit length Overwrite of a 19.00 kfcl square wave 


Density of overwriti are wave (kfcl 


¢ . Left: NLTS dependence on density; Right - Equivalent HTS dependence on density, LF pattern is fixed at 19 kfci, 
overwritten by HF square waves at density 38-160 kfci. These curves demonstrate that the onset of the 
amplification of HTS is similar to the onset of the NLTS in the HF pattern 
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Overwrite: Model of Transition Shifts 


OVERWRITE RATIO = 2 NLTS=30% 


EQUIVALENT TRANSITION SHIFT, NS 


15 20 25 
TRANSITION NUMBER 


Transition Shifts in a HF square wave 
¢ Upper: OW at ratio 2, NLTS=30% 
(© Middle: OW at ratio 10, NLTS =30% 
¢. Bottom: OW at ratio 10, NLTS = 0% 
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Amnplitfication of HTS: Model 


AMPLIFICATION OF HTS: MODELING RESULT 
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¢ Modeling result: Amplification of equivalent HTS with density 
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Overwrite and NLTS: summary 


¢ Overwrite depends both on the LF pattern density (proximity effect) 
and the HF pattern density (NLTS in the overwriting pattern) 
¢ NLTS in the overwriting pattern cause degradation of overwrite 
¢ Overwrite degradation can be explained by oscillations of transition 
shift in the HF pattern. These oscillations are caused by interactions 
between HTS and NLTS 


° Degradation of overwrite 1s nes at high overwrite frequency ratios 
Up to 6-8 dB) 
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Shape Distortions and Random Noise 


¢ PRML Assumptions: 


e Shape of Isolated Transition is Pre-Determined 
e Magnetic Recording Channel is Linear 


— If these assumptions are met, performance of PRML channel is determined only by 
random noise (medium and electronics) 


— If these assumptions are not met, performance of PRML channel is determined also 
by: : 


= Shape Distortions (Equalization) 
-¢ Channel Non-Linearity (NLTS, Partial Erasure) 
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Error Rate and Noise 


[= On-track | | SAM PLOT 
2- misequalized : : 
3 - off track 


ERROR RATE 


25% 50% 45% 100% 
MARGIN, % 


m Noise decreases the slope of the SAM plot 
m Shape distortions: shift of SAM plot 
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Shape vs Random Noise: Distances in the Trellis 


~ 4000 
800 


600 


ERROR RATE 


400 


33% 66% 
-1 O +4 MARGIN % 


¢ Only noise distortions: peak at “0”; convolve noise distribution with this peak 


¢ Shape Distortions: Extra peaks are generated. Each Extra peak is convolved with the noise 
distribution. SAM shifts to the nght 
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NLTS and Partial Erasure 


_ NLTS and PE decreases amplitude of dipulse 


These distortion match the most dangerous error pattern for the ML 
detector. E.g. typical PR4 error sequence is {1,0,-1} 
NLTS 1s “worst case” distortion 
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NLTS Precompensation 


- Left - histograms before precompensation, center - histograms after precompensation, right - SAM 
plot before and after precompensation | 


¢ NLTS is small (about 15% of the bit period) 
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~ NLTS Precompensation 


‘eae Sequenced Amplitude Margin 


e 1 - Random pattern, all adjacent transitions are precompensated the same way 


¢ 2- dibit pattern, same precompensation as (1) 


¢ 3-Random pattern, 2-nd,3-d and 4-th adjacent transitions have different values of precompensation 
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NLTS and Hard Transition Shift 


41039.3 
| 4012.2 
36935.3 
39611.0 
2831.4 
35209.8 
BINS 
3 30808.6 
4623.6 
264074 
» ws 72006.1 
1641574 17604 9 
(BIL ea 
8207.9 $802.5 
41039 | 401. 
0. | —_ 0.0 
“f° 0 +| | n +1 


¢ Left: Pattern is Precompensated, AC-erased medium 


¢ Right - same precompensation, after medium is DC-erased.. NLTS interacts with hard transition shift 
and samples are distorted a | 
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e 1 - Random Pattern, no precompensation 


e 2-random pattern, AC-erased medium and 
precompensation 


e 3-Same as (2) after medium is DC-erased 


| °  4- Medium is DC-erased, each transition is 
i precompensated depending on its polarity. 
Hard and Hasy Transitions are 
| precompensated ‘separately 
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PRML and Non-Linear Distortions 


¢ Error Rate is degraded 
° Precompensation may not work at all 


es At high recording density interactions between NLTS and Hard 
Transition shift are significant . 


Partial Erasure can not be precompensated 
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Pattern Dependent Equalization 


1 - Isolated pulses, optimal equalization ne Sequenced Amplitude Margin 


2 - equalization is same as (1), random pattern 
3 - Equalizer:is adjusted for random pattern 
4 -.same equalization’as (3), isolated pulses 


Which Equalization is optimal? 
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“e alized Method for Measuring Read-Back Non-linearity Using a Spin-Stand 


A new method for measuring the read-back transfer func- 
tion of a magnetic recording system is proposed. For this 
method square wave patterns are recorded over a range of 
linear densities. The resulting read signals should exhibit 
certain superposition properties if the read process were lin- 
ear. The transfer function of a MR head distorts read-back 
signals and makes the read process non-linear. The inverse 
transfer function is modeled using an arbitrary function such 


as a polynomial or cubic spline. This model is then used to. 


linearize the distorted read-back signals. Superposition prop- 
erties are used to evaluate the accuracy of this correction and 
optimize the parameters of the model. The transfer function 
is then obtained directly from the optimum inverse transfer 
function. The effectiveness of the proposed method is con- 
firmed by comparing the read-back transfer function obtained 
experimentally using this method and a MR read head with a 
direct measurement of the transfer function of the same read 
head. 


LT. INTRODUCTION 


The response of an MR head is intrinsically nonlinear. 
By careful biasing and the use of media with a reduced 
Mré product, the severity of this non-linearity can be 
reduced. In practice, significant non-linearity often re 
mains. Partial response maximum likelihood (PRML) 
detection schemes are able to tolerate a large degree of 
inter-svmbol interference by assuming linear superposi- 
tion of signals. Non-linearities in the write process such 
as non-linear transition shift (NLTS) and hard transition 
shift (HTS) seriously degrade the performance of PRML 
channels and these effects have been studied by several 
authors [1,2]. Recent studies have shown that the non- 
linear transfer function of the MR read element can also 
degrade PRML performance (3,4]. 

In this paper a new method for measuring the MR 
transfer function is presented. This method can be im- 
plemented in-situ in a drive or using a spin-stand. This 
method is a generalization of an existing method [5] 
which measures departures from linear superposition of 
Square Wave patterns in the frequency domain. 


| “IL. METHOD ain? 


In order to ee the behavior of the read process, non- 
linearities in the write process must first be controlled. 
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When a square-wave pattern is recorded, NLTS affects 
each transition equally and thus is not observed (6]. In 
this work, hard transition shift is minimized by erasing 
the medium with a very high frequency pattern prior to 
writing. Partial erasure and transition broadening are 
minimized by writing only at moderate densities. For 
linear read-back the signal is then given by the convolu- 
tion sum (after [7]): 


v(t) = )_ afijui.(t — iT) (1) 
where u(t) is the ideal, linear read-back signal, u;,(t) is 
the read-back pulse from an isolated transition, a[t] is the 
sequence of written data taking values of +1,0,or—1, and 
T is the bit period. Consider the following square-wave 
data sequences 


~ afi} = (..+1, 0, 0,-1, 0, 0,+1, 0...}, (2) 
al See le le ele (3) 
= a +a,{i + 2] + afi - 21, oS 4. (4) 


Thus a high frequency pattern of alternating transitions 
can be expressed as the sum of three lower: frequency 
patterns of alternating transitions. From eqs. 1 and 4 


v3(t) = UV) (t) + vy (t i 2T ) + vy (t ae 2T) (5) 


When an inductive read head is used, playback is linear 
and eq. 5 holds as shown in figure 1. However, when a 
MR head is used, this superposition property no longer 
holds as shown in figure 2. 
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3 
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3 

So 

0.5 aes } | | | 
e 100 200° ~=—SSSsétimee (ns) 
V(t) —— v(t) —= vv, (tv (t+2T)+V,(-2T) - - 


me - 


BIG: 1.-Read Signals obtained axperitaéneally using an in- 
ductive head (averaged over.1000 sweeps). “Read-back with ar 
inductive head is a linear process and superposition applies 
Note that the high frequency signal is ae to the's sum 0! 
three low frequency signals. 
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0 100 200 
v(t) —— 94) — 


time (ns) 
V(t)+¥,(t+2T)+¥,(t-2T) - — 


FIG. 2. Read signals obtained using the test MR head (av- 
eraged over 1000 sweeps). 


If we assume the transfer function of a MR head is a 
memoryless function f(z) we can write 


u(t) = f(u(t)) (6) 


where 0(t) is the distorted read-back signal. If we assume 
further that f(z) is invertible then equation 5 can be 
expressed in terms of the distorted signals 


III. EXPERIMENT 


A 3.5”, 1 memu/cm? disk was used with a highly 
asymmetric MR head to test this method. All mea- 
surements were made using a Guzik 1701MP spin-stand 
and 1601+PRML read-write analyzer with optional high- 
speed analog to digital converter. 

High frequency square wave patterns were recorded 
near the inner diameter (I.D.) with densities of 50, 57.5, 
60, 62.5 and 65 kfci. Corresponding low frequency square 
wave patterns were written with densities between 16.6 
and 21.6 kfci. Recording near the I.D. ensures that the 
preamp bandwidth is much greater than the bandwidth 
of the read signals. Gathering data across a range of 
densities makes this method more robust, reducing sen- 
sitivity to noise and residual write non-linearities. Read- 
back signals were digitized at 2Gsamples/sec and aver- 
aged over 100 sweeps. Each erase-write-read cycle was 
repeated 10 times, to minimize the effects of instanta- 
neous flying-height variation, and the resulting signals 
aligned and averaged. The test process is automated 

d all measurements can be made in approximately 2 


f7"(03(t)) = f7" (bi (t)) + f(t + 2T)) + fT ‘(0 (t — 2 nates A polynomial with four degrees of freedom was 


(7) 


this relationship is the key property we use to find the 
MR transfer function. We can model the inverse trans- 
fer function using some general function such as a poly- 
nomial or cubic spline. The problem’ then reduces to 
one of model -fitting- where the model parameters are 
optimized so that eq. 7 is most nearly satisfied. For 
this method we add additional constraints by assuming 
that there is no output from the head for zero applied 
field and the head has unity gain for small signals thus 
f(0) = 0,df/dz|z,=9 = 1. For example a very simple 
model, with zero value and unity slope at the origin is 


“l(r) sa2+azr? + Br? (8) 
Now a this into eq. 7 yields 
U3(t) + avd3(t) + Bd; (t) 
= 64(t) + a0? (t) + Bo (t) 
+0, (t + 2T) + ad? (t + 2T) + Bo3(t + 2T) 


+0, (t — 2T) + ad? (t — 2T) + Bo3(t — 2T) +. (9) 


where € represents the model fitting error. Applying 
eq. 9 repeatedly for different values of t yields an over- 
determined systém of linear equations. This system is 
then posed in matrix form and solved by pseudo-inverse 
to find the ate of. a and B which Minimize. the mean 
squared error < €? >: 

This method is very sensitive to errors in DC levels. 
For ideal linear read-back the .signals.are. symmetric and 
thus the median value of the signal is zero. An invertible, 
memoryless function does not change the rank. order of 
samples and thus the median value of the distorted signal 
is simply f(0) = 0. This property allows the DC levels 
of signals to be set accurately. 


used to model the inverse transfer function. A system of 
equations was constructed from the data and solved to 
yield the four parameters of this model. Figure 3 shows 


the signals of figure 2 after applying the inverse transfer 
function. 


0 100 200 
y(t) ——— v,(t) — 


time (ns) 
V(t)+v, (t+2T)+v, (t-2T) = 


FIG. 3. MR read signals after applying the inverse MR 
transfer function. Note that the high frequency linearized 
signal is equal to the sum of three low frequency linearized 
signals and that peak amplitudes are symmetrical. 


The MR transfer function can be obtained by plot- 
ting the inverse transfer function and exchanging axes 
as shown in figure 4. MR head output as a function 
of applied field can be measured directly by applying a 
known field to the head in a solenoid. There is reason- 
able agreement between the transfer curve obtained us- 
ing this method and the transfer curve measured directly. 
Discrepancies between. the two. methods are greatest at 
around 50 Oe where there appears to be significant hys- 
teresis-and an abrupt change in slope. Hysteresis vio- 
lates our assumption of a memoryless non-linearity and 
a smooth function such as a polynomial can not faithfully 
reproduce an abrupt change in slope. 
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FIG. 4. The MR transfer function. The dashed lines were 
obtained directly by applying a known field to the head with a 
MR head tester. The solid line was obtained using the method 
described here using a 4th order polynomial to model the in- 
verse transfer function. The field strength above the medium 
was inferred by comparing peak signal amplitude with the 
direct measurements of the MR transfer curve. 


The choice of function used to model the inverse trans- 
fer function should have only a small effect on MR trans- 
fer curve obtained by this method. Figure 5 shows MR 
transfer curve obtained for the same system using differ- 
ent functions to model the inverse transfer function. For 
this combination of head and medium it appears that 
choice of polynomial order has relatively little effect. A 
cubic spline, being continuous only up to the first deriva- 
tive, seems to follow the abrupt change in slope more 
faithfully than polynomial models. 
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FIG. 5. MR_ transfer: curves ‘ebtained using dideecht func- 
tions to model the inverse transfer function.:. 

Our A perenes with this method suggests that models 
with more degrees of, freedom, such as this cubic spline, 


may perform better than simple polynomials for good 
data but that where there is significant measurement 
noise or residual write non-linearity a simple model is 
more reliable. 


IV. CONCLUSION 

In this work we present a new method for measuring 
the read-back transfer function of a magnetic recording 
system. In particular this method can be used to find 
the shape of the MR transfer function. Trials with a 
test head have demonstrated results consistent with the 
MR transfer curve obtained by directly applying a known 
field to the head. 

This method may provide a useful tool for in-drive lin- 
earization of read-back signals for improved channel per- 
formance. This method can also be used to eliminate 
the effect of read-back non-linearity in standard tests for 
write non-linearities such as NLTS and partial erasure. 
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Abstract—Many standard tests for write nonlinear- 
ities such as nonlinear transition shift and partial era- 
sure assume a linear read process. Nonlinear effects 
from readback with a magnetoresistive head can make 
the results of these tests difficult to interpret. We 
show how a polynomial can be used to approximate 
the inverse transfer function and thus linearize the 
read channel. Write nonlinearities can be measured 
more easily by performing tests on the linearized chan- 
nel. We also show how linearizing the channel for 
PRML detection is limited by amplification of noise 
over parts of the inverse transfer curve with high dif- 
ferential gain. 


I. INTRODUCTION 


Magnetoresistive (MR) read heads are rapidly replacing 
inductive heads in high-density digital magnetic recording 
applications. The response of MR heads to applied field 
is intrinsically nonlinear. By careful design, good fabri- 
cation and through the use of media with a reduced Mré 
product these nonlinear effects can be minimized. In prac- 
tice significant nonlinearities often remain. 

Many standard tests for write nonlinearities assume a 
linear read process. These tests attribute all observed 
nonlinearity to effects such as nonlinear transition shift 
(NLTS) and partial erasure. Because of this these tests of- 
ten produce unexpected results when used with MR heads 
(1],{2]. In this paper we show how the response of a MR 
head can be linearized and how this procedure greatly sim- 
plifies the interpretation of some standard tests for NLTS 
and partial erasure. We briefly examine PRML detection 
with a linearized channel. 


Il. METHOD 


We assume that the transfer function of the MR head is 
memoryless, continuous and invertible. Under these con- 
ditions the inverse transfer function can be approximated 
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by a polynomial. In this paper a polynomial of the form 
fUl(z) =2+an? + Bx? + yx‘ + 62° (1) 


was used. The parameters a, §, y and 6 of this polynomial 
model were chosen to most nearly linearize the read chan- 
nel. The method for finding these polynomial coefficients 
is described elsewhere (3]. Fig. 1 shows the approximate 
inverse transfer function obtained and the corresponding 
approximate forward transfer function. 

Read signals were linearized by digitizing the waveforms 
at 2 Gs/s and applying the polynomial of eq. 1. The re- 
sulting linearized waveforms were analyzed numerically. 
Time domain methods such as dipulse extraction [4] and 
correlation [5] were performed directly on the sampled 
waveforms. Frequency domain methods such as harmonic 
elimination tests for NLTS [6] and third harmonic ratio 
tests for partial erasure [1] were implemented using a dis- 
crete Fourier transform to extract the desired frequency 
components from the signal. 


Approximate Transtler funcuion 
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Fig. 1. Polynomial approximation for the inverse transfer function. 
Axes are scaled so that the height of a positive isolated peak is one 
unit. 


Ill. EXPERIMENT 


A 3.5", 1 memu/cm? disk was used with a  soft- 
adjacent-layer (SAL) biased MR read head. Read current 
was set at the maximum recommended value of 12mA 
to yield maximum signal amplitude. At this bias cur- 
rent peak asymmetry of approximately 30% was observed. 
Measurements were made using a Guzik 1701MP spin 
stand and 1601+PRML read-write analyzer with optional 
high speed analog-to-digital converter. 


Uncorrected Linearized 
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Fig. 5. Read signal autocorrelation function for a pseudorandom 
pattern. Nonlinear readback produces a strong “echo” which over- 
laps the echo due to NLTS. Linearizing the channel eliminates this 
“echo”. Density is 110 kfci. 


recorded, with precompensation, at a density of 100 kfci. 
The read signal was digitized at 2 Gs/s, linearized, and the 
data clock was recovered by correlating the received signal 
against the known data sequence. The waveform was dec- 
imated using the recovered clock and the resulting wave- 
form was equalized to a PR4 target using a minimum- 
mean-squared-error 24 tap filter. For comparison the 
same proceedure was repeated without linearization. 


Optimal equalization for the linearized signals achieved 
a smaller r.m.s. error and narrower distributions as seen 
in Fig. 6. Fig. 7 shows the sequence amplitude margin 
plot obtained using the method of error filters [9]. The 
SAM error rate for the linearized signal is significantly 
better than for the uncorrected signal for margin up to 
40%, reflecting improved equalization after linearization. 
At higher margin, detection using the uncorrected signal 
performs better than detection using the linearized signal. 
Some parts of the inverse transfer function have a steep 
slope to compensate for the flat response of the MR head 
as it begins to saturate. By a first order approximation 


L f(z) +n 
g{Z) x +ng'(£) 


2 


where x is the ideal undistorted read signal, ¢ is the noisy 
and distorted read signal, f(-) is the MR transfer func- 
tion, g(-) = f~'(-) is the inverse transfer function and n 
is additive noise. Thus steep regions of the inverse trans- 
fer curve have high differential gain which significantly 
amplifies noise. This amplified noise degrades error per- 
formance at high margin. While MR distortion degrades 
the performance of PRML channels it is not possible to 
fully recover this loss by linearizing the read channel. 
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Fig. 6. Histogram of sample values for recording at 100kfci. Lin- 
earization before equalization produces a slimmer distribution. 
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Fig. 7. Sequence amplitude margin at 100kfci with and without lin- 
earization. Linearization before equalization shows improved per- 
formance for small margins. 
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A. Measuring Partial Erasure 


Partial erasure can be measured by recording a square 
wave pattern at density f and again at 3f [1]. Partial 
erasure reduces the amplitude of the high density pattern 
and for linear readback the amplitude loss factor a due to 
partial erasure is given by 
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At low densities, high fields above the medium begin to 
saturate the MR element and thus reduce the amplitude of 
the low density pattern[2]. This effect causes the value of 
a to become negative, incorrectly suggesting a negative 
partial erasure. At very low density both low and high 
density patterns are attenuated by read nonlinearity and 
Qa remains at unity. By linearizing the read channel the 
true partial erasure can be measured as shown in Fig. 2. 
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Fig. 2. Partial erasure amplitude loss factor measured using the 
third harmonic ratio method. Negative values of amplitude loss 
factor at moderate density are an artifact of nonlinear readback. 
Linearizing the channel eliminates this effect. 


B. Measuring NLTS by harmonic elimination 


NLTS can be estimated by writing a pattern which has 
no energy at some harmonic frequency kf. For a linear 
channei the read signal will also contain no signal at kfo. 
NLTS can be inferred from signal amplitude at kfo un- 
der the assumption that NLTS is the only nonlinearity 
present {i], [6]. Nonlinear readback will also contribute to 
the signal at k fo, making it difficult to estimate NLTS in 
this way. Fig. 3 shows net NLTS as a function of write 
precompensation. Ideally net NLTS should approach zero 
for optimal precompensation. When a MR head is used, 
NLTS measured in this way remains high regardless of 
precompensation. Once the read channel is linearized the 
net NLTS drops close to zero for optimal precompensa- 
tion. 
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C. Measuring NLTS by Dipulse Extraction 


By calculating the dipulse response of the channel us- 
ing a pseudorandom test pattern it is possible to measure 
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Fig. 3. NLTS estimated by 5th harmonic elimination. For the un- 
corrected channel net NLTS appears to remain high regardless of 
precompensation. Bit period is 36 ns. Density is 90 kfci 


many nonlinear effects in the channel [4]. Nonlinearities 
in the read process produce several additional “echos” in 
the dipulse response. Some of these additional “echos” 
overlap “echos” for NLTS and hard transition shift and 
make it difficult to interpret the dipulse plot as shown 
in Fig. 4. Linearizing the read channel eliminates these 
additional “echos” . 
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Fig. 4. Dipulse response of the uncorrected channel shows many 
additional “echos” due to nonlinear readback. Density is 110 kfc. 


D. Measuring NLTS by Correlation 


The channel autocorrelation function for a pseudoran- 
dom input is the autocorrelation of the dipulse response 
calculated from the same pseudorandom input. It is thus 
possible to identify many of the same “echos” seen in the 
dipulse response. These “echos” in the autocorrelation 
function can also be used to measure write nonlineari- 
ties [5]. The presence of read nonlinearities produces ad- 
ditional “echos” which overlap “echos” produced by write 
nonlinearities as seen in Fig. 5. Again linearizing the read 
channel eliminates these additional “echos”. 


E. Linearizing the channel for PRML detection 


Nonlinear readback has been shown to degrade 
the performance of partial-response maximum-likelihood 
(PRML) channels [7], [8]. In this section we examine the 
performance of PRML detection on the linearized chan- 
nel. A pseudorandom data pattern of length 127 was 


